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We study theoretically the linear conductance of a quantum dot connected to ferromagnetic leads. 
The dot level is split due to a non-collinear magnetic field or intrinsic magnetization. The system 
is studied in the non-interacting approximation, where an exact solution is given, and, furthermore, 
with Coulomb correlations in the weak tunneling limit. For the non- interacting case, we find an 
anti-resonance for a particular direction of the applied field, non-collinear to the parallel magneti- 
zation directions of the leads. The anti-resonance is destroyed by the correlations, giving rise to an 
interaction induced enhancement of the conductance. The angular dependence of the conductance 
is thus distinctly different for the interacting and non-interacting cases when the magnetizations of 
the leads are parallel. However, for anti-parallel lead magnetizations the interactions do not alter 
the angle dependence significantly. 



I. INTRODUCTION 

Spin transport through hybrid structures of both 
semiconductor and metal systems is now a well- 
established field of researchi. In recent years, also 
spin transport through quantum dots formed by nan- 
otubes, molecular systems or nanoparticles has started 
to attract attention both experimentalljiSiSiiiSiSiLS and 
theoreticallyi^i^°i"i^^i"i^'^i^^i^^ One of the motivations for 
this has been the possibility of using quantum dots with 
a single spin as qubits, but naturally also a number of 
fundamental questions arise on how the interplay of spin 
coherence and interactions influences the transport. 

When the spin coherence is maintained during the elec- 
tron passage through the mesoscopic structure, the trans- 
port needs to be described in a coherent fashion at least 
for spin degree of freedom. Even for weak tunneling one 
therefore has to modify the usual rate equation approach 
if the magnetization directions are non-collinear ^tiii^ 
For stronger tunneling coupling also many of the non- 
equilibrium Green's function methods fail to describe the 
case when both interactions and spin coherence are im- 
portant, such as f or exam ple the Hubbard I approxima- 
tions used in R,efs. ll5llfiL a 

For the strong coupling particular interesting 

situation occurs when the mesoscopic system has a mag- 
netization which is non-collinear with respect to that of 
the ferromagnetic leads, which is what is studied in this 
paper. We study the conductance of a single-level quan- 
tum dot connected to ferromagnetic leads, and, in addi- 
tion, the spin states of the quantum dot being split due 
to, e.g., an applied magnetic field. Experimentally, one 
could realize this geometry if the magnetic leads are fab- 
ricated as magnetic thin films, where the magnetic field 
is strongly pinned in the plane parallel to the film, so 
that a small magnetic field component can be applied 
to the quantum dot without changing the magnetization 
direction of the leads. The geometry of the device is il- 
lustrated in Fig. 1. We describe the system using a mod- 
ified Anderson Hamiltonian where the leads are assumed 
to be polarized and the magnetic field gives rise to an 




FIG. 1: Schematic drawing of the quantum dot system, con- 
nected to two ferromagnetic leads. In addition, there is an 
applied magnetic which spin polarizes the quantum dot in 
a direction non-collinear to the magnetization of the leads. 
In this paper, we focus on how the conductance depends on 
the angle (j>, and the interplay between interactions and spin 
interference. 



energy splitting of the dot level. This gives an effective 
double-slit geometry with two different paths through the 
central region. In the case of parallel magnetizations of 
the leads and no Coulomb repulsion on the dot, we find 
sharp anti-resonances in the conductance for certain an- 
gles (j), which is explained as destructive interference. If 
the Coulomb repulsion on the dot is increased, a cross- 
over to a simple spin-valve behavior is seen. Changing 
the angle in an experimental setup can therefore indi- 
cate whether interactions on the dot are important or if 
a single-electron picture is sufficient. On the other hand, 
in the case of anti-parallel magnetizations of the leads 
no change in the angular dependence is observed when 
the Coulomb repulsion on the dot is increased. Finally, 
we limit ourselves in this paper to the study of linear 
response, so non-equilibrium effects associated with spin 
accumulation are not included. 

The paper is organized as follows. In Sectional we set 
up the model Hamiltonian. We start by calculating the 
conductance for the non-interacting case in Section IIIII 
and in Section llVI we consider the cotunncling limit. Fi- 
nally, in Section we discuss the result and its limita- 
tions. 
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II. MODEL HAMILTONIAN 

The model Hamiltonian of the quantum dot coupled 
to magnetic leads is 



H = Hrr, + HT+HD, 



(1) 



where 



a— L,_R,fca" 

where a denotes the left or right electrode and where we 
have assumed that the polarizations of the electrodes are 
parallel or anti-parallel. The magnetization of the quan- 
tum dot is not necessarily parallel to those of the leads 
and in the spin basis of the leads, the dot Hamiltonian is 



where is the orbital quantum dot energy, B represents 
the magnetic field splitting, and U is the Coulomb en- 
ergy for double occupancy. In a diagonal basis, the dot 
Hamiltonian is 



(4) 



where the d and d operators are related by the unitary 
rotation 



a' 

Finally, the tunneling Hamiltonian is 



cos((/)/2) sin((/)/2) 
-sin(0/2) cos(^/2) 



(5) 



a—L,R kaa' 

Here we allow for the tunneling matrix element to be 
spin dependent, because the states in the leads depend 
on spin direction. The possibility of spin flip processes 
during the tunneling process is not included explicitly 
here, however, it could easily be included by replacing 
the diagonal matrix by a non-diagonal tunneling ma- 
trix. Depending on parameters this would correspond to 
having an angle between the lead magnetizations, which 
would modify the details but not the general behaviors 
that we discuss. 



III. RESONANT TUNNELING WITHOUT 
CORRELATIONS 

We start by discussing the current through the system 
in absence of correlations, i.e. J7 = 0. This treatment is 
relevant where mean-field approaches are valid (of course 



^0 should be determined self-consistently) , e.g. for a reso- 
nant tunneling junction system. The conductance follows 
from the Landauer formula as 



G=- / dtT(e) 



dnp{e] 



(7) 



where T(e) is the total transmission coefficient for both 
spin directions, and np is the Fermi-Dirac distribution 
function. In terms of the retarded and advanced Green's 
function the transmission coefhcient isii 

r(e)-Tr[G''(6)r«(e)G'-(e)r^(e)] , (8) 

with (using the spin basis where Hd is diagonal) 

G^^^ie) = (G„-i - S2'^(e) - S^"(e))"' , (9a) 

[Go(e)],,, = 6a,a'ie -Co- aB)-\ (9b) 

ka" 

where g'2%A^) = (e - C,fc^ ± and r„ = 

—i — SJ^]. In order to simplify the analysis, we as- 
sume that parts corresponding to the real parts of g"^'"" 
give constant shift of and hence can be incorporated 
in Go, and, furthermore, that the tunnel width functions 
^L,R are constant in energy. With these assumptions 



r (A ^L,R f 1 + Pl,r cos (j) Pi.fl Sin 

1 L,r[£) ^ ' 



where 



2 V ^L,flSin0 l-P^„cos( 



(10) 



kcr 



and Pa denotes the polarization of the tunneling from 
lead a defined through ^ = 5 (1 + crPj'^l- Finally, 
we have 



(12) 



With these formulae it is straightforward to find the 
transmission coefficient. 

In the following, we study the current using the sim- 
plifying case of fully polarized leads, i.e. half-metallic 
contacts. For the parallel configuration, Pl = Pr = 1, 
we obtain for the transmission function 



Tp{e) 



^ R 



[x — COS ( 



(1 — x^)^ + {x — cos (j))-^ ' 



(13) 



where 



e - ^0 r" + F" F" , , 
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FIG. 2: (Colors online.) The conductance for the non- 
interacting case in both the parallel (left panel) and the anti- 
parallel (right panel) cases for different values of the magnetic 
strength and at the symmetric point = 0. For the parallel 
case, destructive interference leads to an anti-resonance for 
= 7r/2. The anti-resonance sharpens as B/P" decreases. 
The conductance is maximal when the applied field is paral- 
lel to magnetizations in the leads (<^ = or tt). In contrast, 
for the anti-parallel case the current is maximal for = 7r/2 
and vanishes at = and tt because of the spin- valve effect. 



Let us look at the symmetric condition x — 0, where 
the incoming energy lies symmetrically between the spin 
split levels: 



^ R 



1 



cos^ I 



(15) 



This yields a sharp anti-resonance when the dot magneti- 
zation is perpendicular to the magnetization direction of 
the leads. The anti-resonance occurs due to a destructive 
interference between the transmissions through the two 
spin split levels, as explained in the following section. It is 
interesting to note that the zero transmission at cf) ~ tt/2 
is not broadened by higher order tunneling events, as one 
might have expected. Moreover, for the non-symmetric 
situation x 0, the anti-resonance remains, but is shifted 
away from (f> = 7r/2 to the point where x = cos0, which 
follows from Eq. (|13l) . The anti-resonance found here is 
in fact similar to an effect exploited in single atom trans- 
mission through optical lattices.^® 

Next, we look at the antiparallel situation. Therefore 
we set Pl = ~Pr = 1 and find 



TApie) 



4r|r|sin2, 



D 



(16) 



where 



= 4(1- x^f + 2vlV^ - 



[Vl - Vr? cos^ (17) 



with ya = r'^/B. Note that for the symmetric coupling 
r° = r^j, the only angle dependence is the sin^ (j) factor 
in the numerator. In Fig.|21we show examples of the angle 
dependence of the conductance for both the parallel and 
antiparallel cases. 

A difference between the parallel and antiparallel ge- 
ometry is the existence of an optimal i3-value in the 
antiparallel case. Considering the symmetric condition 



0, we find for the parallel case {4> = 0, tt) 



T'lnax 
P 



(Co) 



pO pO 
^ R 



+ (r72)2' 



(18) 



which decreases when increasing B, because the levels are 
moved away from the chemical potentials of the leads. 
For the antiparallel geometry, the maximal transmission 
at the symmetric condition is (</> = 7r/2, 37r/2) 



AP (4oj 



pO pO 
^ R 



[B+Tirymr 



(19) 



and the optimal value -Bopt = V^^l^^J"^ found. This 
value is a comprise between the spin blockade for B = 
and the levels being far from the chemical potentials of 
the leads for S/T" > 1. 

Finally, Fig. 13 shows how the conductance changes 
when ^0 is varied. In the parallel case, the destructive 
interference only occurs close to the symmetric condi- 
tion, X w 0, whereas when both levels are above or below 
the resonance a simple spin valve effect is observed. If 
one of the levels is exactly at the chemical potentials of 
the leads the conductance is constant (^q = ±B). For 
the antiparallel case no qualitative change in the angular 
dependence is seen. 



IV. COTUNNELING CURRENT FOR THE 
INTERACTING CASE 

In the presence of interactions the interesting interfer- 
ence effect discussed above can be significantly changed. 
Naturally, this is in general a difficult problem to tackle. 
Here we look at the case where the position of the levels is 
far from resonance as defined below, and a second order 
perturbation calculation suffices. This is the "cotunnel- 
ing" regime, where an electron is transferred through the 
system in a two-electron process. The cotunneling app- 



proach is valid far from resonance, i.e. when \e 



Ml 



and + U\ are larger than the bias, temperature and 
tunneling broadening. 

For the resonant case the cotunneling result diverges 
and there is a cross-over to the sequential tunneling 
regime, see App. ^ In these points interactions are not 
important, because transport is determined by only one 
level. 

We calculate the current using a scattering formalism 
and the transition rate per unit time between an initial 
state i (with energy Ei) and a final state / (with energy 
Ef) is calculated as 



27r 



rf, = -\{f\m\'SiEf-E,), 



where the transition operator T is 



T = Ht 



1 



E,-Ha + 111 



-T 



(20) 



(21) 
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FIG. 3: Contour plots of the conductance in the non- 
interacting case for the parallel (top panel) and antiparallel 
(lower panel) geometry for different values of the angle (j> and 
the orbital dot quantum dot energy ^o/ro {B /Tq = 1 in both 
figures). For the antiparallel geometry the same angular de- 
pendence is observed for all ^oi whereas a change is seen for 
the parallel geometry as explained in the text. 

and Hq — Hj^jj + The current is now given by 

processes which bring an electron from the left to the 
right lead, F^j,^, minus the opposite processes: 

J = ei^nL - ^lr)- (22) 

The transition rates will only be calculated to lowest non- 
vanishing order in the coupling. Thus, in the initial state 
the leads can be considered as two non-interacting elec- 
tron gasses (here assumed fully spin-polarized) with dif- 
ferent chemical potentials. The probability to find the 
lead a in a state is denoted W^^ . Because the leads 
are non-interacting, W^^ satisfies 

Three different regimes are considered (see Fig. 0J): 
i) the level with the lowest energy is always occupied 
and the other energy level always empty, i.e. e-^ < ii and 



^2 




(A) (B) 



FIG. 5: The two possible second order processes contributing 
to the tunneling amplitude, here shown for the parallel config- 
uration in the first regime. Because of interactions, the path 
to the right is suppressed when the lower energy spin-state 
is occupied by one electron. Hence the interactions tend to 
destroy the interference effects seen in Fig. |21 valid for U = 0, 
and instead one gets the angular dependence shown in Fig. |S] 

ii) both levels are empty, i.e. £| > /i and > /i, 
Hi) both levels occupied, i.e. < /i and + U < fj,. 
We show how the result is obtained for the first regime 
with parallel geometry and only state the results for the 
other regimes. 

In case i), the initial states can be written as 
I*) = ki-^^fl-T), with energy = E^^ + £'^^ + 
and probability W^^W^^. First, we outline the deriva- 
tion for F^^ in the parallel configuration. In the fully 
polarized case, the leads contain only electrons with 
one kind of spin, say spin-up, and the final states are 
\fkk') = j,|C^ |z), meaning that an electron is 
moved from the left to the right lead in the process. In- 
serting the expressions for the final states and the tun- 
neling Hamiltonian in Eqs. (|20|l and H21|l . we see that 
the first non-vanishing term is second order in the cou- 
pling Hrp. Moreover, two processes give a contribution: 

A) first the electron on the dot jumps to the right lead 
and later an electron from the left lead refills the dot, or 

B) a left lead electron enters the dot, creating a double 
occupied state, and later leaves the dot to the right. The 
two processes are sketched in Fig. |S| The scattering rate 

can now be written as 
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where 



(25) 



(26) 




x(/, 



The matrix elements and are easily calculated, 
and after inserting the coupling constants in Eq. and 
assuming them to be independent of energy, we arrive at 

1 r 1 r 



RL 



R 

27r 



X np(e-A*L)[l - '^pie-Mi?,)]- 



(27) 



For F]^]^ exactly the same calculations are carried out 
and wc finally obtain for the current 



j: 



de 



R |2 



;i 27r 



(28) 

At small bias voltage, we use that jej ^ |e||, + C/| and 
when inserting the form of R from Eq. jsj, we obtain 



J' 



cos2(0/2) sin2(0/2) 



Co-i? ■ ?o + 5 + C^J ' ^^^^ 
where /i^ — /i^ = eV . Note that the two terms in the 
square brackets are the amplitudes for the two tunneling 
processes shown in Fig. 13 In Eq. (|29(l it is evident that a 
destructive interference between the two paths can occur 
at certain angles. 

Of course, a similar calculations can be performed for 
the antiparallel configuration, and with the same con- 
straints the result is 

n 2 



J 



('0 

AP 



p2yiLLR 

Ah 



1 



1 



B + U 



sm^<j). (30) 



Looking at this expression, it is clear that the construc- 
tive interference for the anti-parallel case is not sensitive 
to the angle, because both paths have the angular depen- 
dence. 

The cotunneling results for regime i) are plotted in 
Fig-Elfor different values of the Coulomb repulsion U . It 
is easily verified that to second order in F and U = Q the 
cotunneling results (|29|) and (|30|l and the exact results 
for [/ = derived above coincide. 

In the other two regimes similar calculations give 

"cos2(0/2) sin2(0/2)^' 



J 



Co 



B 



Co 



,■^ F^F" 



Co + S Co - 



B 



sm 



(31) 
(32) 



FIG. 6: (Colors online) The conductance for the interact- 
ing case in both the parallel (left panel) and the anti-parallel 
(right panel) cases for different values of the Coulomb inter- 
action and for the symmetric point = 0. For the parallel 
case, we see how the anti-resonance at (/> = 7r/2 is destroyed by 
interactions, due to the blocking of one of the paths in Fig. |^ 
For large U one instead has a cross-over to a spin- valve effect 
at = TT, where the dot is occupied by an electron with spin 
opposite to the lead electrons. The angle dependence thus 
dramatically depends on the interaction. This is not so for 
the anti-parallel case shown in the right panel, where merely 
an overall suppressing is seen. 




FIG. 7: (Colors online.) The cotunneling result for the paral- 
lel geometry in the three different regimes [j ~ i, ii, Hi] with 
i) 5o = 0, ii) Co/-B = 2 and in) ^o/B = -4. 



for both levels being empty, and 



h 

/ N F"f0 



H 2 



sin (h sin (h 



^o + B + U ^a-B + U 



, (33) 



, (34) 



for both levels being occupied. Note that for [/ = 
holds J*-*^ = j'"-' = and it agrees with the non- 

interacting results in Sec. Illll to lowest order in F^^F^j. 

For the antiparallel geometry the same angular depen- 

(i) -2 

dence is found for all three regimes with G'i'p cx sm 0. 
Fig. 13 shows the conductance for the parallel geometry 
for the three different regimes. It is seen how the destruc- 
tive interference only occurs in regime i), i.e. when the 
levels are split on each side of the chemical potentials and 
the interactions are not too strong. The angular depen- 
dence for the two other regimes depends on which level 
is closest to the chemical potential of the leads. 
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V. DISCUSSION AND SUMMARY 

For the non-interacting case we found the exact hnear 
response conductance. For the parallel configuration we 
saw that the angular dependence could be understood 
in terms of an interference effect occurring due to the 
spin-splitting of the dot. For = the spin bases of 
the dot and the leads are identical, and the electrons can 
only pass through the dot via the spin-t level. For a 
non-coUinear magnetic field the dot spin eigenstates are 
superpositions of the lead spin states, which gives two 
possible paths through the dot. This leads to destruc- 
tive interference at (f) = tt/2, and the resonance sharp- 
ens as the tunneling coupling increases. This result is 
in agreement with the cotunneling result in Eq. H29|l set- 
ting U = 0. Interestingly, the destructive interference 
is destroyed by interactions, giving rise to an interaction 
induced enhancement of the conductance. Thus the con- 
ductance goes from being proportional to cos^ (f) in the 
non-interacting case to being proportional to cos^(0/2) 
for the strongly interacting case. However, if the leads 
are not fully polarized the angular dependence persists, 
but is gradually washed out for decreasing polarization. 

For the antiparallel geometry the situation is very dif- 
ferent. In the non- interacting case (and symmetric cou- 
pling) the conductance is proportional to sin^ for all 
values of B/T^. This is simply understood as a spin- 
valve effect, giving zero conductance at = and (j) — n, 
whereas for non-coUinear i?-field the spin can fiip when it 
passes the dot. Furthermore, in this case there is a con- 
structive interference at </> = tt/2. Now, with interactions 
on the dot one spin channel becomes partially blocked, 
but the angular dependence remains the same even in 
the limit of large U. This is thus in stark contrast to 
the above mentioned cross-over observed for the parallel 
configuration. 

We have thus focused on two cases: 1) strong tun- 
neling and no interaction, and 2) weak tunneling with 
interaction. The cross-over between these two regimes is 
an interesting and challenging issue, because a formalism 
that captures both coherence and correlations on equal 
footing is needed. 



APPENDIX A: SEQUENTIAL TUNNELING 
LIMIT FOR PARALLEL MAGNETIZATION 

In the main part, we have discussed the effects of higher 
order tunneling where interference effects take place. For 



completeness, we here give the results for the sequential 
tunneling limit for parallel magnetization. 

In the parallel case, we can use the conductance for- 
mula derived by Meir and Wingreenii for proportional 
coupling if we also assume that P = Pj^ — Pj^. With 
these assumptions, we have 



h 



where 



|PjTr[r(G»-G^H)]. 

(Al) 



^1 + r°« 



1 + P cos (j) P sin (f) 
P sin A 1 — P cos ( 



To lowest order in the tunneling, which is the so- 
called sequential tunneling case, the Green's functions 
in Eq. (jAl|) are given by 



z(G«-G^)=^'^i'^T' + (l-"i)'5T 











where 



and 



(A2) 



1 + e-/5(2«o-h(7) +2cosh{l3B)e-K» 



(A3) 



(A4) 



where /3 is the inverse temperature. Thus for the parallel 
case, the conductance becomes a sum of four peaks: 

ro^ro = {h{^, + B)[l-n,]+h{^, + B + U)n{) 
X (1 + Pcos<j)) 

+ (/i(Co -B)[l~ n^] + h{^o -B + U)n^) 
x(l-Pcos0), (A5) 

where h{uj) = — -^j-^^ • This agrees with the non- 
interacting result for J7 = and F" 0. 

For the antiparallel geometry it is more complicated 
to obtain a simple result because the current formula in 
this case contains the lesser Green's functioniii The full 
rate equation results have been given in Ref.m 
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